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Abstract 

Particularities of the statistical mixed phase model for describing the deconfinement 
phase transition and physics constraints on the model parameters, are discussed. The modi- 
fications proposed concern an improved treatment of hadron-hadron interactions within the 
nonlinear mean-field model of nuclear matter and an inclusion of the one-gluon exchange 
correction into the quark-gluon sector of the mixed phase thermodynamic potential. It is 
shown that the mixed phase model successfully reproduces both the available lattice data 
on the QCD deconfinement transition and thermodynamic characteristics of hadron systems 
evaluated within the modern models. The 'softest point' of the equation of state resulting in 
a possible formation of a long-lived fireball in nuclear collisions (the 'softest point' effect) is 
analyzed. The obtained results are confronted with the data on the two-phase model. Some 
manifestation of the mixed phase formation and deconfinement transition is estimated in 
fireball expansion dynamics. We discuss experimental possibilities of observing the softest 
point effect as a signal of the deconfinement phase transition. 
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1 Inroduction 



By now, a search for signals of the QCD phase transition from a hadronic state to that of 
unbound quarks and gluons (quark-gluon plasma) has a long story. As was noted many years 
ago, the presence of a phase transition in the equation of state (EoS) might essentially influence 
the evolution of a highly heated and compressed system (a fireball) formed in relativistic heavy-ion 
collisions. In particular, in the vicinity of the phase transition point EoS is getting softer, which 
results in a qualitative change in the hydrodynamic picture of the fireball expansion. Interest 
in this has been recently revived in respect of extensive discussion ||, |^, H, ||] of the so-called 
'softest point' effect (SP effect) which is expected to effect a general strategy of the experimental 
search for QCD phase transition signals. 

The SP existence is related to a local minimum in the pressure-to-energy density ratio p/e as a. 
function of energy density e and gives rise to elongation of the expansion time of the hot and dense 
nuclear matter created or, in other words, to formation of a long-lived fireball [|I], ^. Qualitatively, 
it is understandable in the following way. If the fireball formed in heavy-ion collisions is in a state 
close to the thermodynamic equilibrium, the subsequent stage of collective expansion needs time 
to be inversely proportional to the sound velocity squared. For the equilibrium nuclear matter 
the last quantity is close to the p/e ratio and coincides with it in the ultrarelativistic limit. 
Hence, when the energy density inside the fireball reaches the value at the SP Ssp its life time is 
getting longer due to the local minimum in p/e at e = Ssp- 

The most detailed and reliable lattice QCD calculations are available only for the case of 
baryonless matter, = 0. They show that a sharp increase in the energy accompanied by a 
relatively small increase in pressure occurs within a quite narrow interval of temperature near the 
deconfinement point T^ec @, H; H- Therefore, one should expect that experimental manifestation 
of the long-lived fireball formation may be observed in a rather narrow range of colliding energy 
Eiab around the value which corresponds to reaching esp- However, a reliable estimate of this value 
of Eiab corresponding to the SP effect is not simple due to two types of available uncertainties. 
The first type is related to the choice of EoS resulting, in particular, in different values of Ssp] 
the second is caused by uncertainties of the hydrodynamical description of the fireball expansion 
(a different degree of energy thermalization in the initial state, inclusion of dissipative effects and 
so on). In this paper, we would like to present a statistical mixed phase model for deconfinement 
transition and demonstrate the reliability of the proposed EoS for the thermodynamic description 
of nuclear matter and the SP effect as a peculiarity of this EoS. We also look in a simplified 
qualitative way at the dynamic manifestation of this peculiarity during the nuclear fireball decay. 

As a rule, EoS with the QCD phase transition is treated phenomenologically within a two- 
phase model which assumes that the system studied may be realized in either hadronic phase or 
quark-gluon plasma 0]. In such a consideration for the case of baryon-rich matter, the concavity 
condition of the free energ with respect to volumy is broken in the region of phase transition. 
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This leads to the necessity of introducing the Gibbs mixed phase where spatially separated quark- 
gluon plasma and hadrons coexist The presence of the Gibbs mixed phase smooths a sharp 
change of thermodynamic quantities in the deconfinement region but only for ub 7^ 0. As to the 
case of non-zero baryon density, there is a jump by an order of magnitude in average energy 
and entropy at the transition point T^ec- However, as has been recently shown by the lattice 
calculation for the SU{3) system with two- flavor quarks |^, the QCD phase transition at ra^ = 
is a continuous crossover with gradually increasing energy and entropy densities in the region of 
0.8 Tdec <1.2 Tdec [ii S S Hi- should emphasize that this result concerns only the ub = 
case. At present, there is no fully satisfying generalization of the lattice approach to the case of 

7^ 0, just because the phenomenological statistical models are quite appealing to get unique 
information on the deconfinement transition in baryon-rich matter. To get rid of discrepancy 
with the lattice results, some attempts are made to use various smoothing procedures for the 
energy (or entropy) jump obtained in the two-phase model [|I], 0, ^, the smoothing region being 
called a generalized mixed phase in order to distinguish it from the Gibbs mixed phase to be 
realized in the framework of the two-phase model for ub = only at Tdec- Arbitrariness of the 
smoothing procedure changing from paper to paper as well as its artificial character result in a 
large uncertainty of estimating esp which varies from 0.7 GeV/ fm? to l.'bGeV/ fm? H, |^, ^. 

To avoid the above mentioned difficulties, our study of the EoS features related to SP is based 
on the statistical mixed phase model developed in [|1^ rather than on the two-phase approach. 
In our model, the role of 'the generalized mixed phase' to smooth sharp transitions in the system 
is played by a mixed state of nuclear matter where unbound quarks and gluons coexist with 
hadron phase without any spartial 'stratification' or bubble/drop formation. In other words, 
average distance between hadrons and quarks/gluons may be of the same order as that between 
hadrons. In contrast with the Gibbs mixed phase, in our case the hadron interaction with unbound 
quarks/gluons plays an important role and should be taken into account, which makes the problem 
be more complicated. Below, this state of nuclear matter will be referred to simply as a mixed 
phase (MP). As has been shown in Ref. |]12|, in the region of QCD phase transition a MP state 
turns out to be thermodynamically more preferable as compared to both the pure hadron state 
and quark-gluon plasma, all sharp changes in thermodynamic quantities to be specific for the two- 
phase picture of deconfinement being smoothed noticeably in the presence of MP. In particular, 
for the = case of S't/(3) symmetry with quarks of two light fiavors the crossover transition 
was obtained at the (pseudu) critical temperature Tdec = 150 MeV [0], i.e. one succeeded in 
getting agreement not only in the value of deconfinement temperature Tdec consistent with lattice 
predictions in the interval T^ec = 143 -r- 154 MeV but also in the type of phase transition. 

The paper is organized as follows. In Sect. 2, the structure of the statistical model of MP is 
outlined and its modifications and improvements made after the first publications are discussed in 

*This system is a good approximation for the real case. 
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more detail. New lattice data for thermodynamics of the SU{3) system with two light quarks are 
confronted with calculational results within the MP model. Specific features of SP in the mixed 
phase are studied in Sect. 4: the value of esp, the depth of a local minimum, the baryon-density 
dependence of esp and so on. Expansion dynamics for a fireball described by the MP model with 
deconfinement transition is investigated in Sect. 5. In the concluding part, the results obtained 
are summarized and their possible relation to experiment is discussed. 

2 Mixed Phase Model 
2.1 Model structure 

At small energy density (lower than that inside a nucleus in the ground state), a nucleon is a good 
quasiparticle defining thermodynamic properties of nuclear matter. When the energy density 
increases, other hadron species should be involved into consideration, first pions and delta^. 
Hadron abundance is determined by minimization of the free energy of the system with respect 
to the hadron number what gives rise to certain relations between chemical potentials for different 
species of quasiparticles named chemical equilibrium conditions. This allows one to present 
all thermodynamic quantities as a function of only two variables: Temperature T and baryon 
density ub (or baryochemical potential fi) which are just the most widely used thermodynamic 
variables fixing a macroscopic state of nuclear matter^. At an essential increase in the energy 
density, nuclear matter is expected to suffer a phase transition into quark-gluon plasma where 
new quasiparticles, unbound quarks and gluons, are now playing a decisive role. Thermodynamic 
properties of a fully equilibrium quark-gluon plasma depend on T and as well. Till the recent 
time the two-phase model, as has been noted in the introduction, remains as the main tool of 
studying the transition from a hadron phase to a state of unbound quarks and gluons. In this 
approach, two completely different thermodynamic potentials are defined: One is for hadron phase 
fhadiT) and the other is for the quark-gluon plasma fqcpiT). At the given baryon density, one 
may find a temperature T{nB) at which both the potentials are equal to each other. The curve 
Tins) is a borderline for this phase transition. Above this curve, the quark-gluon plasma turns 
out to be preferable since its free energy is smaller than that for a hadron phase. Respectively, at 
temperatures below T^np) a hadron state is preferable. It is quite evident that the structure of 
the two-phase model predetermines the first order of the QCD phase transition. Indeed, since the 
ways of determining the functions fhadiT) and fqcpiT) are quite different and independent, it is 
unlikely that the equality of free energies at the transition point fhad(T{np)) = fqcpiTinp)) will 
be accompanied by the equality of their derivatives f'had{T{nB)) = fqcpiTinB)), saying nothing 

^We have in mind here quasiparticles with in-medium dispersion relation rather than free hadrons 
*In this paper we are limited ourselves by consideration only a system with two light flavors. 
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of coincidence of higher derivatives. Therefore, a typical situation for the two-phase model is 
characterized by the equality of thermodynamic potentials but not their derivatives. It is just the 
case of the first order phase transition. On the other hand, it is clear now why the application of 
the two-phase model to deconfinement results in the lack of agreement with the lattice calculation 
to predict a continuous phase transition for the SU{3) system with two light flavors. 

Instead minimizing the two different free energies separately at the given values of T and 
and then choosing the one that is smaller at the chemical equilibrium point, it seems to be more 
efficient to construct a unique thermodynamic potential depending on a number of both hadrons 
and unbound quarks/gluons. The minimization of this potential with respect to all quasipar- 
ticles will determine their concentration. In other words, in the treatment of thermodynamic 
behavior of highly excited nuclear matter one should allow quarks and hadrons to coexist in a 
phase-homogeneous way and in a chemical equilibrium state. The question of whether such a 
physical state is possible or not is solved by studying the stability of a mixed phase. As has been 



shown in [|T2[ and mentioned above, in the deconfinement region MP is thermo dynamically more 
preferable than both the hadron phase and the quark-gluon plasma. Thus, taking account of MP 
is inherent in the consistent treatment of the QCD phase transition. In addition, the scheme 
considered does not predetermine the order of a phase transition. 

To find the free energy of the mixed phase we shall use the following effective Hamiltonian ||12|| : 

H = I ^+(5,s)(v^+f/,({p}))^,(f,s)cix 

+ [ {^-V^ + ml+U,i{p}))^P,ix,s)dx+ (1) 

q s •' 

+ ^ni(^' ^) (\/-V2 + + Unj{{p]) ) ^„j(x, s)dx- 

nj s 

- C{{p})V, 

where ipa(yX,s) denotes a field operator for quasiparticles of a— species: a = g means unbound 
gluons; a = g is unbound quarks of the q -type { q = u, u, d, d, ■■■ a = nj correspond to 
n— quark (gluon) hadrons of the j-type. Index s accounts for quantum degrees of freedom (spin, 
isospin, color). In (|l]) Ua{{p}) is the mean field by which the system acts on a quasiparticle of the 
a— species with the mass nia and {p} is a joint set of individual densities pa of all quasiparticles. 
It is worthy to note that C{{p}) V is the c— number term, where V is the system volume. The 
separation of quasiparticles implies the passage from an exact Hamiltonian depending on field 
operators of 'generic' particles ip to an approximate one to be written in terms of the quasiparticle 
operators ip- 

As a rule, the quasiparticle vacuum, satisfying the condition ^|0 >= 0, does not coincide with 
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the particle vacuum to have iplO >= . In other words, 



< 0|/7|0 >7^< 0|/J|0 >~< 0|iJ|0 > . 

By the vacuum definition < 0|if|0 >= and therefore < 0|-?/|0 >^ 0. It is easy to see that the 
last relation would be broken if the correcting c— number term was absent in the Hamiltonian ([^). 
Now the physical meaning of this term is getting clear: This is the system energy (taken with the 
minus sign) in the case when its microscopic state is the quasiparticle vacuum. 

One should clarify why [/^({p}) and C{{p}) do not depend explicitly on temperature. The 
point is that one should be very careful in working with the Hamiltonian to be dependent on 
thermodynamic variables like quasiparticle density. If one succeeds in separating quasiparticles 
accurately starting from the exact Hamiltonian, no problem with the thermodynamic consistency 
arises. Such an example is given by the Bogoliubov model of the weakly interacting Bose gas [|13|. 
However, this can be done not in every case, and then the effective quasiparticle Hamiltonian 
should be constructed under some assumptions. In this case, it is quite probable to meet the 



problem mentioned. As has been shown in Refs. |12, 14|, Hamiltonians depending on thermody- 
namic variables must satisfy certain relations. In the opposite case, one may meet non-equivalence 
of different statistical ensembles and finally thermodynamic inconsistency when different calcu- 
lational ways of some thermodynamic characteristics lead to different results. Relations, which 
should be taken into account while constructing an effective Hamiltonian, may be called conditions 



of thermodynamic consistency and written down as follows [0, |T^ : 



(_)=0. {^J=0(Va), (2) 

where (A) denotes the average value of the operator A over statistical ensemble. With the 
Hamiltonian (|1]) conditions (0) may be presented in the form |n 



^ dUa dC ^ dUa dC 

It is easy to convince ourselves that eqs.(^ are compatible only if the mean fields do not 
depend explicitly on temperature. Really, the second equality in may be rewritten as 

J2 PaUa - C = (p, 
a 

where (j) = 4>{{p}) denotes some arbitrary function of quasiparticle density. Then, by differenti- 
ating this expression with respect to pp 

^ — UfS -\- 2^ Pa ^ - ^ , 

dpf5 a dpp dpfs 
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and using the first equation in d), we get 

g0({p}) 



/3 ! 



which proves the temperature independence of Up. One should note that this consideration does 
not exclude completely an explicit temperature dependence of any mean field. We are dealing 
only with the semi-relativistic Hamiltonian (|l]) where the energy of free quasiparticles is given 
in a relativistic form but their interaction with surrounding matter is described by means of a 
non-relativistic mean field. We shall come back to the discussion of this approximation later on 
when some specific parametrizations of the mean fields are considered. One should note that the 
conditions of thermodynamic consistency for the case of the relativistic scalar field were discussed 
in Ref. It is of interest that here an opposite situation is realized: The mean field must be 
explicitly temperature dependent. 

The thermodynamic consistency conditions (|^) lead to one more interesting relation to be 
used below. To get it, let us differentiate the first equation in (j^) with respect to the density ps 

dpf3 a " dpsdpp dpsdpfs 

Interchanging the indices 6 and f3 and comparing the obtained expression with the previous one, 
one finds 

dUs _ dUp 



dpp dps 

to be valid if their mixed derivatives are equal, i.e. 



(4) 



dpsdpfs dppdps dp&dpp dppdps 

This is fulfilled if the second derivatives of Ua and C are continuous, which takes place almost in 
all cases, in particular, in the mean field approximations used below. 

So, we have discussed the structure of the effective Hamiltonian for the mixed phase. Now 
one should specify the mean field acting on a quasiparticle. 

2.2 Interaction of unbound quarks/gluons with the mixed phase 

The mean field affecting color charges in the quark-gluon plasma is used to be taken inversely 



proportional to some power 7 0.3 -i- 1 |16, O] of the density of surrounding color particles 



The simplest arguments in favor of such an approximation are the following. According to the 



MIT- model of hadrons |]T8[, the average energy of a system of unbound quarks and gluons may 
be presented as 

Epi = E Pg(^)^ + EE sqrte + ml + BV , 
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where Pg{k) and Pq{k) are the momentum distributions of quarks and gluons and B is the so-called 



vacuum pressure. As has been proposed in ||T^, this expression may be rewritten as follows 

Here ppi denotes the density of all partons in the quark-gluon plasma 

ppi = Pg + J2 PQ- (5) 

So, it is reasonable to take the single particle spectra of quarks and gluons in plasma in the 
following form: 

ujg{k) = k + —, ujq{k) = Jk^ + ml + — (6) 
Ppi Ppi 

where the parameters A and 7 may be found, for example, by fitting the available lattice results 
on QCD thermodynamics(see ||12|| ). 

At the same result for u)g{k) and u!q{k) one may arrive from a completely different consideration 
based on string potential arguments |jT6|. An essential point is an assumption on the nearest- 
neighbor interaction. It allows one to avoid divergences since in this case the integration of the 
string potential y ~ err is carried out over a sphere with the radius to be equal to the average 
distance between the nearest neighbors rather than over the whole space. There are two arguments 
in favor of why it is possible to limit ourselves to accounting only the interaction between the 
nearest neighbors. One of them (see the first two cited papers in |jl6|) is that the screening of a 



point-like color charge has an analogy with the Debye-Hueckel and Thomas-Fermi screening in 
an ordinary plasma. A common feature here is a correlation between charge density fluctuations 
in different space points: The correlation is essential at a short distance between these points but 
it is negligible at a large one. Thus, the interaction of two plasma constituents situated far from 
each other will be proportional to the product of their average color charges which are, however, 
equal to zero because the quark-gluon plasma is color neutral, locally and on the whole. The 
second argument is based on the consideration of a probability p{r) to form a string of the length 
r (see the last cited paper in |TB[) which turns out to be an exponentially decreasing function. 
As a result, we have that at the calculation of the interaction energy of a quark with surrounding 
plasma one should integrate not the string potential but an effective one V ~ p{r) ar to fall 
down quickly when the distance exceeds the average distance between the nearest neighbors. 

One should notice a specific feature of eqs. (|]) to be directly related to the confinement 
phenomenon. The last implies that an infinite energy should be spent for appearance of an 
isolated quark or gluon. The case of isolated color charge corresponds to the limit Ppi 
and, as is seen from (^, the energies ujg{k) and ujq{k) are really going to infinity in this limit. 
Therefore, singularity in spectra (P) at the zero density Ppi ensures confinement of color objects. 
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On the other hand, deconfinement possibihty is provided by the fact that ujg{k) and ujq{k) are 
reduced to non-interacting particle spectra at high values of the energy density (or ppi). 

It is important to explain yet one more peculiarity of (^. The point is that the terms A/p^i 
and ruq + A/p^i are frequently interpreted as effective masses of unbound gluons and quarks. 
But under such an interpretation we have an unusual dispersion relation. These constructions 
are named semi-relativistic ones because they combine the allowance for a non-relativistic 
interaction with a relativistic form of the kinetic term. Their implementation is argued by the fact 
that the term A/p^i noticeably exceeds the average momentum of unbound quarks and gluons 
in the region of thermodynamics variables where the QCD phase transition takes place {ppi is 
small). Consequently, the application of (|) is expected to be equivalent here to making use of 
the following spectra: 



UJg{k) 



\ 




Ppi 

The interpretation of A/pJ^i and niq + A/ppi as effective mass of unbound gluons and quarks could 
be quite satisfactory if a new version of dispersion relations, ujg{k) and ujq{k), would allow one to 
develop a thermodynamically consistent approach. However, the statistical model based on such 
spectra can not be thermodynamically consistent (see Appendix A). The solution of this problem 
is that the mean field in the new version of dispersion relations should be a function of 'scalar' 
densities rather than 'usual' particle densities, similarly to the relativistic scalar self-consistent 
field in the Walecka model |2^. And so, the relativistic analog of (^ should be taken as follows: 



Ug{k) 



k^ + 



Ps,pl 



Uq{k) 



\ 



k"^ + \mq + 



A 

7 

Ps,pl , 



(7) 



where 



Ps,pl = p. 



s,9 



s,q 



and, by analogy with the Walecka model, the quantities ps^g and ps^q may be called scalar densities 
of unbound gluons and quarks 



p. 



A. 

27r2 



A/Pl 



pi 



's,pl 



exp 



A;2 dk 



ujg{k) 



- 1 



27r2 



mq + A/p^ 



s,pl 



P + ( 



T^q + ^/plpl 



exp 



P dk 



UJqjk) - Pq 

T 



+ 1 



Here and pa denote a number of quantum degrees of freedom and a chemical potential for 
quasiparticles of the a— species [pg = 0). No problem on thermodynamic consistency arises at 
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making use of (|^. Moreover, when the average momentum of unbound quarks and gluons is 
small, as compared to their effective masses, the scalar and ordinary densities are close to each 
other, ps^pi ~ ppl , and, therefore, the results of the QCD phase transition consideration based 
on (|) and (|^ should be close too. Thereby, in the first approximation, the quantities A/p^; and 
niq + A/ppi may be interpreted as effective masses of unbound gluons and quarks. 

Up to here, a pure quark-gluon plasma without any hadron admixture was discussed. How 
spectra of unbound quarks and gluons will be look like in the mixed phase (MP) where the 
presence of hadrons influences on the self-consistent field affecting color particles by environment? 
Based on the results for the pure quark-gluon plasma, it is reasonable to consider two limiting 
cases. The first one is to use eqs. for MP as well, which implies neglecting any interaction 
between hadrons and unbound quarks and gluons. It corresponds to so strong boundness of 
hadron constituents that the presence of free color charges in their surrounding does not result 
in their color polarization, i.e. hadrons remain to be color neutral and due to that are unable 'to 
see' quarks and gluons not involved into hadrons. The second limit is based on the relations 

Ugik) =k + ^, uj,{k) = ^k^ + ml + ^, (9) 

where p is the total density of quarks and gluons in MP: 

P = Ppi + Yl ^P^J ■ (^0) 

nj 

This approximation corresponds to a very strong color polarization of hadrons when free color 
charges in the MP make no difference between bound and unbound quarks/gluons. At first sight, 
one can hardly make any choice between these two limiting cases, (^) and (|^), because no estimate 
is available for hadron polarizability in color environment. Fortunately, there exist arguments of 
thermodynamic nature in favor of choosing (^). 

Indeed, a large set of quantum numbers characterizing quark/gluon results in a high variety of 
different gluon and quark-gluon systems suffering the deconfinement phase transition of different 
types. In particular, in the pure gluonic SU{2) system, the deconfinement is the second order 



phase transition with a peak-like behavior of the heat capacity ||2T]|. In the case of the SU{?>) 
system without quarks there is a phase transition of the first order ||22|, while for the baryonless 
SU (3) system with dynamical quarks of two light flavors the deconfinement transition is of a 
crossover type |^, ||, 0. As has been noted in the introduction, there is no satisfying lattice 
results for real baryon-rich > systems yet. So, a natural demand for a statistical model is 
to reproduce the whole variety of known phase transitions. If we keep the relations (^, at the 
zero baryon density we get 



Ppl 



27r2 







fc2 dk 


exp 


\k + A/plA 


- 1 


T 



+ 
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+ y — 



exp 



T 



+ 1 



Equations of this type were investigated numerically (see the second reference in and the third 
one in |T^) and it was shown that their positive solution ppi > exists only at T > and 
appears by a jump. It results in the first order QCD phase transition for both quark-gluon systems 
with an arbitrary number of flavors and a purely gluonic case. In the simplified consideration 
of (|TT|) (classical approximation and rriq = 0) this numerical analysis can be exemplified by the 
equation 

1 A 

Ppi = -i (^9 + E Q e^P( - jr^)- 

^ q '''Ppi 

Under these assumptions the latent heat Ae of the phase transition depends only on the parameter 
7 



23 



£SB \ ^ 37) ( 7 ) 
where £55 is the volume energy density of an ideal quark-gluon gas. The value of Ae may be 
equal to zero only in the limiting case 7 — > but it is not interesting since to have a confinement 
one needs 7 > . Thus, the choice of MP spectra in the form of (^) necessarily leads to the 
first order phase transition independently of the quark-gluon system in question to be in contrast 
with the available lattice QCD data. Therefore, the limiting case (H) corresponding to a strong 
polarization of hadrons in the MP by surrounding color charges is more relevant. As has been 
shown in the earlier study [|12|, the statistical MP model making use of (§) does allow one to 
reproduce a large variety of phase transitions in an agreement with the lattice data. 

Concluding this section, one should emphasize once more two main points essential for getting 
single particle spectra (|^): They are the confinement of color charges and variety of deconfinement 
scenarios. The first results in that the mean field affecting color charges in MP should be taken 
inversely proportional to some power of quark/gluon density; the second needs that it should 
be just the total density of quarks and gluons, both unbound and bound into hadrons. These 
propositions are necessary conditions to get the representation (^ but they are not sufficient 
conditions allowing more complicated approximations for unbound quark/gluon spectra. For 
example, 

A 

ujg{k) = k + 



{ppi + A Y.npnj)^ 

nj 



^ V^+^-^^l^. (12) 

nj 
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where the parameter A ~ 1 specifies the polarizabihty degree of hadrons in MP. Below, a more 
simple version will be used which assumes in the Hamiltonian (|I]) 



(13) 



If a more detailed description is needed, eq.(|l2|) may be used without any trouble. 



2.3 Interaction of hadrons with the mixed phase 

In MP in addition to hadrons there exist unbound quarks/qluons, the average distance between 
the first and last being of the same order as that between neighbor hadrons, as noted above. Due 
to that there is no a well-separated spatial boundary between coexisting phases that is specific 
of the Gibbs mixed phase. Therefore, the proper description of the hadron component of MP 
should take into account not only ordinary and well-studied hadron-hadron interactions but also 
new quark/gluon-hadron interactions. 

The mean field acting on the nj— hadron in MP may be presented as follows: 

f/n, = f/i? + f/S'^ . (14) 

Here, L^^^"* is esual to f/„j in the limit when MP is degenerated with respect to unbound quarks/gluons 
and thereby reduced to a purely hadron state. Such a situation occurs at rather a low energy den- 
sity, beyond the deconfinement region [l^ . In its turn, the term v'^f' is caused by appearance 
of unbound quarks and gluons in the system. In accordance with this definition, we have 

lim f/if = 0. (15) 



It is noteworthy that it would not be completely correct to interpret the first term in (|lj) simply as 
a hadron interaction energy with the MP hadronic component and the second term, accordingly, 
as a mean field acting in the MP on this hadron from the quark-gluon component side. Indeed, the 
interaction nature of hadrons in MP should be changed because of their polarization by unbound 
color charges. Therefore, the second term in (|1^) should include not only the interaction energy of 



a hadron with the quark-gluon component but also some correction to the hadronic term arising 
due to the polarization mentioned. 



It seems to be problematic how to find U'^l^ ■ However, eq. (HI) coming from the conditions of 



thermodynamic consistency (|^) simplifies this procedure. Using (|]) and taking into account that 
uj{f is independent of densities of unbound quarks and gluons, from (|) we get 



dp, dpnj p^+i ' dpg dpnj p^+i ■ ^ ^ 
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The integration of eqs. (^6|) gives 

uif = ^ + ^({p4) (17) 

with an arbitrary function of hadron densities ip{{ph}) (here {ph} denotes a set of hadron densi- 



ties). Then, confronting (KB) with (|T^, one gets 



f4f = ^{l-il-w,,r) , (18) 

where the concentration of quark-gluon plasma in the mixed phase Wpi = Ppi/p may be considered 
as an order parameter of the QCD phase transition in this modeL Thus, if Uq and Ug are known, 
the thermodynamic consistency conditions allow us to find unambiguously the correction term 
(H). It is not surprising because the mean fields acting on unbound quarks and gluons, Ug and Ug, 
carry out information on hadron polarization and their interaction with the plasma component 
of MP to define a form of U^j\ At finding U^jK the conditions of thermodynamic consistency 
allow us to use this information without its explicit separation. 

According to (|T^), the term U^^ is the mean field acting on the raj— hadron in the degenerated 



case of MP, i.e. in the case of purely hadronic environment. In the earlier version of the statistical 



MP model ||Tj, hadron-hadron interactions were treated in the Hartree approximation and 



u!!^^ were taken as follows: 

U^^^ = V $ . ■ n ■ (19) 

mi 

where ^nj,mi is the integral of the interaction potential between the nj—ih and mi—ih hadrons 

^nj,mi J ^nj,mi(^^)dT . 

To get rid of inconveniences to work with a large number of different hadron-hadron potentials, 
scaling relations may be used 

*ni,mi(r) = ^ $„«c(r) . (20) 

Here ^nuc{f) denotes some effective nucleon-nucleon potential. The approximation (|20| ) was ar- 
gued in Refs. |jT2|, |1^ by the consideration of conservation laws in the reaction of fusion of colliding 



clusters like in the molecular dynamics. Similar relations arise when the hadron interaction is 
reduced to the interaction of constituent quarks but in this case it is more preferable instead of 
([20| ) to use an equivalent relation between the hadron coupling constants: 

gnj/gmi ^ n/m . (21) 

It is sufficient to take account of hadron interaction in the Hartree approximation at the zero 
baryon density owing to the high temperature of the QCD phase transition. As to finite baryon 
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densities, the Hartree approximation results, as it is well known, in the lack of important features 
of the thermodynamic behavior related to the saturation point. To get an appropriate description 
of MP in a large range of T and variables, we use in this paper other approach based on the 



Walecka-like |20] mean- field model. At present, there are various modifications of this model. 



Among them, the structure of baryon spectra proposed by Zimanyi et al. [24| is the most suitable 
allowing us as before to work with the Hamiltonian (|T]): 



UJnj{k) = \/k^ + mlj + gr,nj + 9a,nj '^2{y) , (22) 

where x and y are given by 

^ ^ 9r,mi Pmi y V ^ ^ 9a,,mi Pmi ■ 
mi mi 

Here gr,mi > and ga,mi < are the repulsive and attractive coupling constants, respectively, and 
functions ipi{x) and ip2{y) are defined by the equations: 

bi(pi = x, -bi{ip2 + b2^pl) = y (23) 



with two free constants hi and 62- In 01 ^ mixture of nucleons and A isobars was considered 
and the correct value of the ratio for their coupling constants put an additional constrain on the 
conditions ( p2|) and (p3|). In our model, not only baryons but also mesons are included. To reduce 



a number of free parameters it is convenient to use again the scaling approximation (PT|). Indeed, 
using (pTD and (|2^), the hadronic term may be represented as: 

t/g) = n (<^i(p - ppi) + ^2(p - Ppi) ) , (24) 

where ipi and ip2 satisfy the equations 

cii^i=p-ppi, -C2^2- c^^l = P - Ppl (25) 

with p—ppi = J2nj '^Pnj- In (0) Ci, C2 and C3 are the following combinations of the old parameters: 

_hi _hi _ 6162 

Ci — — ; C2 — — , C3 — — — 



9r 9a 9a 

where gr = 9r,nj/n ga = ga,nj/n. Values of the new parameters Ci, C2, C3 are fixed by properties 
of the ground state (T = and = no ~ 0.17 fm^^ = 1.28 x 10^ MeV^) of nuclear matter: 
Pressure is zero, binding energy is —16 MeV and compressibility is 210 MeV. Under these 
conditions we get ci = 5.48 x 10"*^ MeV'^, C2 = 5.87 x 10^ MeV^, C3 = 26.9. 

The advantage of the scaling relation ( pT]) is seen from the compact form of (p^). But it is 



unclear yet what are the limits to applying this approximation. In the hadron sector of MP at 
temperature T < 150 MeV and baryon density ub < 10 uq nucleons, the A isobars and pions 
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dominate in the two-flavor consideration. And so, here one may hmit oneself to the consideration 



of only these three species of hadrons. According to (]2^) , the mean fields acting on nucleons and 
A's from the hadron environment are the same which is not in disagreement with Zimanyi's et 
al. paper where the ratio 1/1.3 was obtained for nucleon and isobar coupling constants. In 
accordance with that, making use of (|2^ ) one may underestimate the interaction energy of isobars 
with surrounding hadrons by about ~ 30% but the A contribution to major thermodynamic 
quantities like pressure or average energy is noticeably smaller than that of pions and nucleons 
even near the boundary of thermodynamic variables T < 150 MeV, ub < lOfiQ. Thus, the 
application of is justified for baryons if bulk thermodynamic properties are studied. One 
should note that it is not the case when the A multiplicity is of interest; then deviations from 
the scaling behavior ( ^Ij) must be taken into account. 

Let us consider now the applicability of the scaling relation (pTf) and, accordingly, the ap- 
proximation (|2^ ) for pions. In the Walecka-type models ||20|, pions are used to be treated as 
an admixture of noninteracting particles to a nucleon liquid. However, simple arguments coming 
from an excluded volume consideration follow a quite different picture: In the dense nucleon 
matter pions strongly interact with the environment. It is confirmed by detailed microscopic 
calculations within the both relativistic virial expansion p5| , |26[| and Bruckner theory that 
turn out to be close to the results obtained within the MP model. Really, let us turn to the recent 



paper by Rapp and Wambach ||2^ where thermodynamics of interacting pion gas was studied in 
the most comprehensive way. In our approach, the hadronic resonances are singled out in advance 
according to a general statistical treatment since Hagedorn |2S|. To compare with Ref. [^, where 
a p— resonance channel of the tttt interaction is taken into account, one should calculate thermo- 
dynamic quantities for a system of interacting pions and p— mesons. This system is governed by 
the same Hamiltonian (|I|) where only the vr— and p— terms should be kept with eq.(P^ used for 
the mean field Unj- With this Hamiltonian the following expression can be derived for quarks 
bound into hadrons: 



P- Ppi 



dk 



exp 



A;2 + mlj + U, 



nj 



T 



(26) 



1 



Note that meson chemical potentials are equal to zero and certainly ppi = 0. If the density p 
satisfying (|2T|) is known, the pressure and average energy for the system of vr and p— mesons may be 
estimated by means of the Hamiltonian (||). At the given mean fields Unj the correcting c-function 
in the Hamiltonian is found by solving the differential equations (§). All formulas to be necessary 
for getting general thermodjTiamics quantities from (|l|) can be specified in Refs. |jl2|, |14[. The 
numerical results obtained with the scaling relation (pT]) for the pressure p energy density e in the 
vr, p— system are presented in Fig.l and 2. For comparison, the results by Rapp- Wambach 



are 



plotted in the same figure for a gas of interacting vr— and p- mesons as well as for non- interacting 
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ones. It is seen that the allowance for meson interactions in any way displays a qualitative 
difference from the case of ideal-gas mesons. It is of interest that at T < 190 MeV the attraction 
plays a more essential role. At these temperatures the average energy is higher than in a mixture 
of non-interacting mesons because the interaction potential is getting negative. Oppositely, at 
T > 190 MeV the repulsion is dominated, which results in slowing down the growth of p and e 
and then they are getting even smaller than the appropriating values for non-interacting meson 



gas. It is seen that the statistical MP model with using (p^ fairly well agrees with the results of 
Ref. at temperatures < 120 MeV. 

Thus, eq. (|2ll) is a good assumption for describing thermodynamic properties of an inter- 
acting meson system. But how is about ttN interaction ? Hamiltonian for the meson-nucleon 
system may be constructed in the same way as in the previous case but nucleon and antinu- 
cleon terms should be added to the vr— and p— ones. The appropriate terms should be included 
also in ^6|). In Fig. 3 the meson potential U2j obtained within the MP model is presented for a 
TT, p, A^, iV— system at = . It is seen that the energy of a pion at rest (i.e. the sum of pion 
mass and U2j) amounts to 108 MeV at T = 150 MeV. For the same quantity at normal nuclear 
density and T = 150 MeV the value 110 MeV were obtained by Shuryak p6| within the rela- 
tivistic virial expansion. Therefore, the approximation (|2^) based on the scaling relation (^) and 



nonlinear Zimanyi's version ||2^ of the mean-field model of hadron interactions well reproduces 
bulk thermodynamic properties of the probe systems considered. 



The relations (0), (pilD , ([T8|) and (p4D describe all the mean fields acting on the MP con- 



stituents and define the volume energy density (with the minus sign) of quasiparticle vacuum C. 
The last quantity is found by solving the set of differential equations (|^) that are convenient to 
be presented as follows: 

dUg ^ dUq> ^ dUnj dC 
dUg ^ dUg ^ dUmi _ dC 

dpnj V ^Pnj mi 9pnj 8 Pnj 

Casting ( p^ ) into (pT]), the first two equations are reduced to 

dC dC . _^ 

dpg dpg 

which allows one to represent the quasiparticle vacuum energy in the form 

C = --^Ap'-^ + M{ph}), (28) 
1-7 
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where ipi is some function of real variables. The substitution of ( pH]) into the last equation of ( pTf ) 
results in the equation 

which may be solved easily as: 

^1 = Y^:^ + {p- ppd^'fiip- ppi) + ^2{p- ppi)) - 

p-ppi 

{(pi{t) + <f2it)) dt + const . (30) 

The constant in (^) should be taken as zero because in the absence of quarks and gluons (both 
unbound and involved into hadrons) the energy of the whole system should be equal to zero. 
Finally, for the energy density of the quasiparticle vacuum we have 

C = A p'-' [l-il- WpiY"'') + 

p-ppi 

+ {p- Ppi){Mp- ppi) + Mp- ppi)) - J {Mt) + Mt))dt. (31) 



The relations (|I3D , (|I4D , (0), (|2^) and (|3ll) complete the definition procedure of the MP 
Hamiltonian. It remains only to point out which quasiparticles are included into our consideration 
and what are the values of the model parameters A and 7. Being limited in the given paper to 
the case of two light quark flavors, unbound u and d quarks together with their antiquarks are 
considered in addition to unbound gluons. The current quark mass rriq is 7 MeV. In the hadronic 
sector, all mesons up to 0(1020) made of u d quarks/antiquarks are taken into account. Besides 
baryons, nucleons and A isobars as well as their antiparticles are included. This set of hadrons 
seems to be sufficient to get necessary accuracy for general thermodynamic quantities in MP. 

As to A and 7 parameters, the structure of the MP model does not leave much space for 
their values. Indeed, the way based on the MIT bag model for estimating Ug and Ug dictates 
that the parameter 7 ~ 1 and should have the same value for both purely gluonic and quark- 
qluon case. Fitting the lattice data on the QCD thermodynamics for SU{2) and SU (3) quarkless 
systems [p!2| , |l^ gives good agreement between the MP model and the lattice results at 7 ~ 
0.6 0.65 for the SU{2) symmetry and at 7 ~ 0.5 -i- 0.65 for the SU{?>) color group, being the 
best at 7 ~ 0.62 in both the cases. Below, the parameter 7 is fixed just by this value. At last, the 
parameter A describing intensity of color interaction is defined by the QCD running constant. In 
the case of massless fiavors, the last is presented as follows 

1 



as{r) 



— Inh^ 
27r VAr 
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where ho = ^ — ^Nf {Nc and Nf denote a number of colors and flavors), A is the known 
QCD parameter. In accordance with this formulae, the transition from the quarkless SU (3) case 
to that of the SU{3) theory with two light flavors is accompanied by an increase in by about 
10%. So, one may expect that quantities depending on intensity of the color charge interaction 
will not change much at a transition like that. A lattice estimate of the ratio \fa jviip [o is the 
string tension, rrip is the p— meson mass) testifies in favor of this conclusion: ^Jo jvTLp = 0.54±0.05 
for SU{3) without quarks, that should be compared with y/a/nip = 0.545 ± 0.026 for the real 
quark-gluon system (see the review article Q). So, it is quite reasonable to use the same value 
of the string tension -y/o" = 420 ± 20 MeV in both the versions of the SU{3) theory, purely 
gluonic and with quarks of two flavors. The MP model analysis of thermodynamics of the SU (3) 
gluonic system showed that ^4^/*^^^+^^ ~ Tdec [|12|. According to the recent lattice results p2|, the 
deconfinement temperature of the quarkless SU (3) system is estimated as Tdec/ = 0.625±0.003, 
which gives rise to T^ec = 262.5 ± 12.5 MeV". Therefore, for the MP model parameters 7 and A 
we have: 



7 



0.62 , = 250 ^ 275 MeV . (32) 



Now the Hamiltonian of the statistical MP model has been defined completely and we may turn 
to the analysis and numerical solution of the basic model equations. 



3 Mixed Phase Thermodynamics 

3.1 Thermodynamic characteristics of MP and their comparison with 
the lattice QCD results 

Thermodynamics of an equilibrium system is completely defined by the thermodynamic function 
depending on thermodynamic variables to be controlled by the appropriate conservation laws. In 
MP in question there is an admixture of unbound quarks and gluons, and one should elucidate 
how color neutrality of the system is controlled. In the general case, it is possible that every 
color-spin projection of an unbound quark of the given species corresponds to its own chemical 
potential which defines the average number of these quasiparticles. At the passage from one 
color-spin state to other, the chemical potential may be changed, then different quark-color states 
are not equally probable and MP is not color neutral. Otherwise, when all color-spin projections 
are specified by the same value of the chemical potential, the unbound quark admixture in MP 
has on average zero spin and no color charge. Therefore, the color neutrality of MP is achieved 
by introducing a unique chemical potential pa for all color a— quasiparticles. 

Let us consider the temperature and energy- density dependence of the average quasiparicle 
density pa- By calculating the statistical average value < ip^i^a > in the standard way for the 



18 



Hamiltonian (|1]), we have 



27r2 



P dk 



exp 



T 



(33) 



Here z/^ takes the value 1 for bosons and —1 for fermions. Eq. (|33|) is apphcable for unbound 
gluons if Trig = is assumed. Quasiparticle densities (|33| ) formally look like independent quantities 
specified by their chemical potential pa- However, a strong interaction between quasiparticles may 
give rise to attainment of not only thermodynamic but also chemical equilibrium, where chemical 
potentials of all quasiparticles satisfy the relation: 



(34) 



Here is the quasiparticle baryonic charge and p is baryochemical potential. As a result, pa is a 
function of only two variables, temperature T and baryon density ub (or baryochemical potential 
p). 

In accordance with expressions for the mean fields (|TB|), (^^, ( [T^ ) and (^) and chemical 
equilibrium conditions (pif), to calculate pa in the equilibrium case by means of eq. (|33|) one 
needs to know the total quark-gluon density p, the density of unbound quarks and gluons ppi 
as well as baryochemical potential p. To define these three unknown quantities there are three 
equations. The first of them is given by (|) relating ppi with the densities Pg and Pg. The second 
is eq. ( |To|) to be the definition of p. Finally, the third expresses the baryon charge conservation 
law: 



(35) 



The numerical study of these equations has showed that for any values of T, ub from the region 
of < T < 500 MeV and < ub/uq < 40 there is the only solution in the whole interval of 
parameter values = 250 ~ 275 MeV . It excludes the possibility of the first order phase 

transition which in the mean field approximation is accompanied by the presence of a branching 
region for the order parameter, to be Wpi = Ppi/p in our case. 

It is of interest to compare thermodynamic characteristics of MP with the QCD lattice results 
for the SU(3) system with quarks of two flavors. For this aim, one should first find the free energy 
of MP. Since the Hamiltonian (Q) is given in the mean field approximation, the free energy is 
calculable exactly 



/ 



E 



2 7r2 



In 



1 — Ua exp 



k'^ + ml + Ua - Pc 
T 



k'^dk + 



+ E ~ ^ ■ 



(36) 
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Then, all other thermodynamic quantities may be easily obtained by the proper differentiation 
of the free energy. In particular, the differentiation of (|36|) with respect to fia results in pa to 
be coincided with expression (P^D- This fact illustrates the thermodynamic consistency of the 
statistical MP model. 

Before proceeding to comparison of our results with the lattice data, one should note the 
following point. At taking into account the interaction of unbound quarks and gluons, we were 
limited by the so-called non-perturbative part of this interaction. There is, however, a pertur- 
bative part corresponding to the one-gluon exchange interaction. It plays no essential role in 
the region of the QCD phase transition and does not change the order of the phase transition. 
Due to that, this term was not earlier included into consideration in Refs. |14[. However, the 
perturbative potential part may noticeably manifest itself at high energy density (temperature) 
and influence some more delicate characteristics like the SP position. In the case of = 0, 
it is possible to take account of the perturbative interaction of quarks and gluons by a simple 
reduction of a number of effective degrees of freedom via the following substitution in the free 
energy density (see PTj]): 



The running coupling constant g may be taken as a constant value because of its weak temperature 
dependence. Fitting the lattice data for the purely gluonic SU{3) theory in the region of 
the QCD phase transition gives g"^ ~ 10. This estimate is expected to be valid also in the SU{3) 
gluon matter with two flavors. So, there are only two free parameters allowed to vary within a 
quite narrow range: 

^1/(37+1) = 250 ^ 275 MeV , g^ ^ 10 . 

As is follows from the analysis of the MP thermodynamics at hb = 0, the best agreement with 
the lattice approach is achieved at = 8 -h 9. The deconfinement temperature is fixed by a heat 
capacity maximum, see Fig.4. The numerical analysis shows that its value in a simple way is 
related to the model parameters: Tdec = A^^^"'~^^ — 97 MeV. The above noted uncertainties in the 
value of results in a temperature spreading of about Tdec = 153-7-172 MeV. It is noticed 

that two schemes, those by Wilson and Kogut-Suskind, usually used in the lattice approach 
for the discretization of the QCD action agree with each other in both the order of the phase 
transition in the SU (3) system with two flavors and the value of the deconfinement temperature: 
The first scheme gives T^ec ~ 150 MeV ||3^; the second one is T^ec = 143 -h 154 MeV (see the 
first reference cited in P]). The best agreement with both the schemes is achieved by the choice 
^1/(37+1) _ 250 MeV for which we have Tdec = 153 MeV. It is the value of the parameter A to 
be used everywhere below. 

One should note that as compared to the previous version of the MP model, this version is 
not only based on an improved treatment of hadron interactions but also uses the new lattice 
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data for the purely gluonic SU{3) case. For the case with quarks, the results of the detailed 
comparison between the MP model and the lattice data are in good general agreement with each 
other as it is seen in Fig. 5 from the presented temperature dependence of the reduced energy 
density and pressure. However, the Kogut-Susskind scheme has a problem with the description 
of hadrons: The energy density and pressure for both the new data obtained on the lattice 
with a number of temporal steps Nt = 6 and the former ones for A^j = 4 |^ go to zero too fast 
as it would be expected for a pion gas at T < Tdec and ub = 0. It is explained (see the forth 
reference in 0) by the fact that for the lattice sizes considered the ratio of the calculated pion 
mass to the temperature of the QCD phase transition is still too high, m^/T^ec = 1.94, though the 
used values of the light quark masses are close to the physical one: rriq/Tdec = 0.075. This small 
paradox is possibly coming from the point that a number of temporal steps in these calculations 
is not large enough. Indeed, the study of a peak of the quark condensate derivative with respect 
to the quark mass in the crossover region shows (see the first cited paper in 0) that changing 
in Nt from 4 to 12 at almost the constant ratio rrig/Taec ~ 0.08 results in decreasing the peak 
height by about twice and extending its width near the ground from 10 MeV to 20 MeV. In other 
words, the crossover is getting to be less sharp when A^^ increases due to a relative increase in 
thermodynamic characteristics below T^ec and, accordingly, to a decrease in the calculated pion 
mass. 

3.2 Softest point effect in the MP model 

The temperature dependence of the energy density e and pressure p in Fig. 5 really presents EoS 
for the SU{3) gluonic matter with quarks of two light flavors at the zero baryon density. These 
results of the MP model together with 'experimental' points of the lattice QCD calculations are 
replotted in Fig. 6 but in the other representation: p/e ratio vs e . One should note that the MP 
results are given here only for the case of = 8, being practically coincident with the curve p/e 
for = 9. The lattice data in the Wilson scheme and the MP model exhibit a minimum in the 
p/e function in the region of the deconfinement transition predicting the softest point position as 
large as esp ~ 0.5 GeV/ fm? . The recent results obtained in the Kogut-Susskind discretization 
scheme plotted in Fig. 5 do not unfortunately allow one to make any definite conclusion as to 
the SP existence and its position due to the above-mentioned problem concerning hadron degrees 
of freedom resulting in underestimation of the pressure at T < Tdec- In Fig. 6 one can see also the 
two-phase model data with the bag constant B = 235^ MeV^ used. This value of the vacuum 
pressure B results in a reasonable value of T^ec ~ 160 MeV but gives a too large value for SP 
esp ~ 1.5 GeV/ fm? to be hardly conciliated with the lattice data predictions. It is of interest 
that an artificial smoothing of the entropy (or energy) jump to signal about the first order phase 
transition within the two-phase model allows one to get results which are very close to those of 
the MP model and lattice QCD if the width spread is about 0.1 Tdec 0- 
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It is noteworthy that the two-phase curve p/e in Fig.6 was used in two independent hydro- 
dynamic calculations and resulted in essentially different estimates of the colliding energy Eiab 
where the SP effect may manifest itself. In the papers by Rischke et al. assuming a full 



stopping, a deep minimum in the excitation function for the directed flow has been predicted at 
Eiah ~ 5 AGeV while a smooth curve corresponds to EoS of a hadron gas. Hung and Shuryak 
assumed that only half an energy is attainable for the hydrodynamic stage that shifts the expected 
energy to a noticeably higher value Eiab ~ 30 AGeV. At this energy and ub = the authors 
obtained almost a double increase in the fireball life-time caused by a very deep minimum oi p/e 
s^t £sp (see Fig.6). In the case of 7^ the two-phase model results follow practically the same 
curve as in Fig.6. But when the baryon density increases, the hadronic states, getting closer and 
closer to the mixed Gibbs state (the hyperbolic part of the p/e curve in Fig.6), are realized keep- 
ing unchangeable the position and general depth of the minimum. One should emphasize that in 
this model there is the first order phase transition and the absence of detailed lattice QCD data 
does not allow one to construct any smoothing procedure to mimic a crossover transition like in 
the Ub = case. However, the statistical MP model developed above may be directly extended 
to the case of baryon-rich matter without any additional parameters. 

As is seen from Fig. 7, the baryon density dependence of the p/e function in the MP model 
differs essentially from the two-phase model predictions: with increasing ub the function minimum 
is washed out and disappears completely at ub > 0.4 no. This leads to quite controversial 
conditions to observe the SP effect in heavy-ion collisions. Indeed, the low value of esp specific 
of the crossover transition testifies that the observation of a long-lived fireball as a signal of the 
deconfinement phase transition might be observed at rather moderate energies, below 10 AGeV. 
However, due to high stopping at these energies, the fireball formed will have a high baryon 
density which, according to the MP model but in contrast with the two-phase model, should 
result in washing out the minimum in p/e and degenerating the SP effect. Some alternative case 
is given by heavy-ion collisions at ultrarelativistic energies of LHC and RHIC colliders. If the 



Bjorken mechanism |^ is realized at these energies, the manifestation of the SP effect is possible 



in the course of evolution of a baryon-free fireball. 



3.3 An estimate of the MP evoluton 



To see experimental consequences of the mixed phase EoS, the latter should be incorporated 
in a hydro dynamic-type model. In the first approximation, we consider an expanding blob of 
the compressed and heated QCD matter (a fireball) formed in heavy-ion collisions. The initial 
state for this fireball is defined based on kinetic results. Calculated within the quark-gluon string 



model the time evolution of ub and e for central Au + Au collisions is shown in Fig.8. The 
expansion stage starts from the time moment when the maximal value of baryon density is reached. 
This instant corresponds to a full overlapping of colliding nuclei. But some later time moment. 
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when the total baryon density in the fireball coincides with the density of participant nucleons^, 
is taken as the initial state for the present calculations. Further fireball evolution is described 
as an isentropic scaled expansion, following Bjorken According to this prescription, the 

volume densities of the conserved quantities, entropy and baryon number, are evolved inversely 
proportionally to the expansion time; so all the other thermodynamic quantities can now be easily 
found for the given EoS at every temporal step (see, for example pT|). 

The results of these dynamical calculations for general thermodynamic characteristics are 
given in Fig.9. The energy density in the fireball is steady decreasing according to the power 
low with the exponent 1.11 instead of 4/3 as expected for the ideal gas case. So, the energy 
density of MP evolves slower. The pressure behaves similarly but it can hardly be approximated 
by a single exponent at Eiab = 150 AGeV. At high colliding energies, the temperature fall- 
off essentially differs from ultrarelativistic ideal gas expansion. The presence of crossover phase 
transition results in a slowdown of the temperature fall-off and is observed as some 'shoulder' 
in the T(t) dependence. One should note that for the first order phase transition the function 
T{t) is not monotonous: the temperature increases in the mixed Gibbs phase This 



irregularity is manifested even in the beam-energy dependence of temperature calculated in the 3D 
hydrodynamic model with EoS which includes the first order phase transition 0]. The proximity 
of the MP evolution to that of the pure hadron phase as well as their difference is understandable 
if one looks at the temporal development of the concentration of unbound quarks and qluons: 
Wpi amounts a noticeable value only in the very beginning of the expansion stage at relativistic 
energies and then decreases exponentially. 

According to Ref. |1[, a pure hadronic phase may be characterized by a certain level of the 
energy density e and this gives some feeling as to the life-time of the fireball formed. In the 
considered energy range the fireball life-time (defined for example by the e = 0.13 GeV/ fm? 
level) increases steadily and exhibits no distinct maximum predicted in [||]. The absence of a 
long-lived fireball is a direct consequence of the mixed phase EoS where the softest point is 
washed out with an increase in ub-, as noted above. It is illustrated again in Fig. 10 where the 
key evolution quantity p/e is shown as a function of the energy density attained in expansion: 
making use of the two-phase model essentially overestimates the SP effect. 

As has been noted, the SMP model predicts that the quark-gluon admixture in the fireball falls 
down almost exponentially with time and its slope strongly depends on the bombarding energy. 
Even near the freeze out point this admixture amounts (1 — 2)% for Eiab > 10 AGeV. On 
the one hand, it means that treating the freeze-out stage within the MP model is not simple and 
deserves special attention. On the other hand, it gives a possible explanation of the recent finding 
that the freeze-out parameters extracted from experiments at Eiab = 10 and 200 AGeV turned 
out to be just on the phase boundary (estimated in the two-phase model) between hadronic and 

^In the kinetic treatment, the participant nucleons are those that have suffered at least one colhsion. 
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quark-gluon phase ||3^ . Since the plasma admixture at Eiab = 3 AGeV is essentially smaller, one 
may expect that the freeze-out parameters will be below the phase boundary in this case. 



4 Conclusions 

So, the hypothesis on coexistence of hadrons and unbound quarks/gluons in nuclear matter has 
been realized within the statistical mixed phase model. The conditions of thermodynamic con- 
sistency as well as general requirements for the confinement of color objects and the variety of 
possible QCD phase transitions constrain the chosen form of the mean fields acting on quasipar- 
ticles in MP and predetermine it, to a certain extent. Making use of the nonlinear self-consistent 
mean-field model to take into account hadron-hadron interactions allows one to describe correctly 
the ground state properties of nuclear matter. Inclusion of the one-gluon exchange corrections 
into interaction of unbound color objects results in the correct asymptotic behavior of thermo- 
dynamic quantities. The developed version of the MP model successfully reproduces the lattice 
QCD results and predicts the crossover-type transition of the deconfinement in the SU (3) system 
with quarks [ub = 0) at T^ec = 153 MeV. Allowing a natural parameterless extension to the case 
of baryon-rich matter, the MP model is in reasonable agreement with the modern-model results 
for describing thermodynamic behavior of hadron systems. 

The MP model predicts that the softest point is located at esp ~ 0.5 GeV/ frn?. This value 
of e is close to the energy density inside a nucleon and thereby reaching this value signals us that 
we would be dealing with a single 'large' hadron consisting of deconfined matter. As has been 
noted above, slower evolution of a system in the SP vicinity as compared to the neighbor regions 
of e may give rise to some observable signs of the QCD deconfinement. Two circumstances 
should be emphasized here. First, the measurements in some definite interval of energies Eiah 
for excitation functions of the quantities of interest are needed, since the crossover transition 
occurs in a finite region of e. Even highly precise measurements at a single energy can not be 
interpreted unambiguously as to the SP effect. Second, in contrast with the two-phase model 
with the first order phase transition, in the MP model the SP effect is washed out at ub > 0.4 no, 
and thermodynamic behavior of the mixed quark-hadron system at T ~ T^^c turns out to be 
closer to thermodynamics of pure hadron phase but not to the two-phase model thermodynamics. 
Therefore, a search for a decisive deconfinement signal implies a comparative analysis of the 
measured excitation functions and some hydrodynamic results based on EoS within both the MP 
and interacting hadron models. 

The directed nucleon fiow characterizing the nuclear matter bounce off in respect of the re- 
action plane is an example of the observable to be used for this kind of analysis |3^. If an 



interacting system lives longer at other equal conditions, it results in loosing the memory on the 
initial reaction plane, i.e. in decreasing the directed particle flow. The most complete analysis 
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of the directed flow excitation function has been carried out in the cited papers by Rischke et 
ah 1^, but with the two-phase model EoS. As has been argued above, the same hydrodynamic 



model but with the MP EoS is expected to give results very close to that obtained in Q for 
pure hadronic EoS. It is of interest to note that though till now there was no direct experimental 
search for the effect in question, the available data at Eiat = 2 and 10 AGeV ||3^ do not appar- 



ently confirm the deep minimum in this excitation function at Ei^b ~ 5 AGeV predicted by the 
two-phase model of the QCD deconfinement. 

The directed flow is certainly not the only signal of the mixed phase and deconfinement 
phase transition. A generalization of the MP model to the case of three flavors shows (see the 
second reference in |Q) a sharp increase in the relative fraction of strange quarks (unbound and 
involved into strange hadrons) to be observed at £ > esp, i.e. the strangeness enhancement is 
correlated with the SP position for ub, this correlation being kept for the case of baryon-rich 
matter. The available systematization of strange particle production may be considered 



as some evidence in favor of such a behavior: the reduced strangeness suffers a sharp increase 
somewhere at Eiab = 5 — 10 AGeV and then practically is independent of the colliding energy. 
Identical particle interference deserves particular attention since it is sensitive to 'granularity' of 
an emission source Alongside with the source radius, the correlation length for unbound 

quarks/gluons will play a role of a new additional scale in the MP case. The production of 
dilepton pairs and hard photons is of most interest. In the course of evolution, a system described 
by different EoS will have different paths in the configuration space, say in (T, fi) space and, 
therefore, different emission probabilities. In addition, taking into account interactions between 
hadron and quarks/gluons will open a new channel for dilepton and soft photon production. 
In particular, this channel may contribute to the enhanced yield of dileptons with the effective 
masses 0.2 — 0.6 GeV recently observed by the CERES Collaboration in relativistic heavy-ion 
collisions Finally, one should note that in the MP model there are no such phenomena as 



overheating and supercooling in the limit of infinite volume if the deconfinement is of a crossover 
type. 

The discussed consequences of a possible realization of MP in heavy-ion collisions are mainly 
of qualitative nature. To get quantitative predictions the developed mixed phase EoS should be 
included into full hydrodynamic calculations. This work is now in progress. 
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Appendix A 



Let us consider a system whose thermodynamics may be descibed in terms of a single type of 
quasiparticles distributed in volume V with density n. Let the Hamiltonin be represented as 
follows: 



H = Y: V^TlPata. - CV, (Al) 



k 
k,s 



where U = U{n) is the density dependent mass of the quasiparticle and the physical meaning 
of C was discussed in the principal part of this paper with respect to eq. (|lD. The creation 
ai'^ and annihilation ag^ operators for a quasiparticle with the momentum k in the quantum 
state s satisfy the Bose or Fermi commutation relations. We want to show that the choice of the 
Hamiltonian in the form ( |A1| ) results in losing the thermodynamic consistency. For simplicity, 
let us limit ourselves to the classical approximation (the high temperature limit). Then, the free 
energy reads as 



f exp(-^^^^^] edk + {fin-C)V , (A2) 



27r2 7 ^ \ T 







where ^ denodes a number of internal degrees of freedom of quasiparticles, /x and T are chemical 
potential and temperature. Using the equality 



2n 

one can get {N = nV): 



n = ^ 7exp (-y^^f±^] edk (A3) 



'dF\ ^ 7 f Vl^TTP\ dU Uk^dk (dc\ 



'dF\ 

Since ( j = /i, the relation ( [A4[ ) is reduced to 

oo 



'8C\ ( fexpf--^) , " "^eik. (A5) 



dn ^ 27r2 7 T / ^fW+lP dn 







Therefore, for the given Hamiltonian C should be a function of both density and temperature 
C = C(n, T) but it leads to some contradiction. Indeed, for the average energy of the system we 
have ^ 

E = (H) = ^ [Vk^ + exp I - ^^'J" ^ - I ^dk - CV. (A6) 
2iT^ J \ T I 







But E may be defined also in another way by means of the well-known formula 

'dF 
df 



dF\ 

E = F-T{ — ) , (A7) 
Ny 
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which gives 



oo 



T 



) 



E 



27r2 



ykP^+lP exp 



dfc - 








n 



(A8) 



Thus, two ways of calculating the average energy of the system, to have the same result in 
both the ways, leads to essentially different expressions, which is considered as thermodynamic 
inconsistency. If the mean field additionally depends on temperature, an extra term proportional 
to the derivative of U with respect to temperature will come into ([A8|). Then, the comparison of 
(|A8|) and ( |A7|) will allow one to get a connection of temperature derivatives of U and C which 
(together with (^)) is a set of differential equations defining such a dependence of on n and T 
under which the conditions of the thermodynamic consistency are fulfilled. 
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Figure captions 



Fig.l Temperature dependence of the energy density of a pion gas (in units of T^) within the 
statistical MP model. Dashed and dotted curves are the Rapp-Wambach results for interact- 
ing and ideal pions, respectively. 

Fig. 2 Temperature dependence of the pion gas pressure. Notation is the same as in Fig.l. 

Fig. 3 The energy of a pion at rest in a pion-nucleon mixture at hb = uq. The curve is calcu- 
lated within the MP model, points are taken from Ref. . 



Fig. 4 Temperature dependence of the reduced heat capacity for the SU{3) system with light 
quarks of two flavors at ub = {cv,sb is the heat capacity of quark-gluon ideal gas). Curves 
1 and 2 are calculated for g"^ = 8 and g"^ = 9, respectively. The deconfinement temperature is 
shown by the arrow. 

Fig. 5 The reduced energy density and pressure of the SU (3) system with two light flavors at 
ub = calculated within the MP model for g"^ = 8 (curve 1) and g^ = 9 (curve 2). Triangles and 
squares are the lattice QCD results obtained in the Wilson and Kogut-Susskind schemes, 
respectively. 

Fig. 6 The p/e-representation of EoS for the two flavor SU{3) system at ub = 0. Triangles are 
the lattice data in the Wilson scheme 0, the solid and dotted lines are results of the MP and 
two-phase models, respectively. 

Fig. 7 Energy and baryon density dependence of the p/e ratio in the statistical MP model. 

Fig. 8. Evolution of baryon number and energy density in central Au + Au collisions at different 
bombarding energies as calculated in the quark-gluon string model. Both quantities are referred 
to the cylindrical {R = 4 fm, L = 2R/'jcm) Lorentz contracted cell (fireball) in the center-of-mass 
frame. 

Fig. 9. Evolution of different thermodjTiamic quantities for the fireball formed in central 
Au + Au collisions at the beam energy Eiab- Calculations are carried out with the MP model EoS. 

Fig. 10. The correlation between p/e and the energy density reached during expansion of the 
fireball formed in central Au + Au collisions at Eiab- The results are obtained within the MP (the 
left-hand side) and two-phase (the right-hand side) 
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